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Abstract 

The temporal evolution of microwave pulses transmitted through random 
dielectric samples is obtained from the Fourier transform of field spectra. 
Large fluctuations are found in the local or single channel delay time, which 
is the first temporal moment of the transmitted pulse at a point in the output 
speckle pattern. Both positive and negative values of local delay time are 
observed. The widest distribution is found at low intensity values near a phase 
singularity in the transmitted speckle pattern. In the limit of long duration, 
narrow-bandwidth incident pulses, the single channel delay time equals the 
spectral derivative of the phase of the transmitted field. Fluctuations of the 
phase of the transmitted field thus reflect the underlying statistics of dynamics 
in mesoscopic systems. 
41.20.Jb, 05.40.+j, 71.55.Jv 
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Statistical optics has concentrated on fluctuations of reflected and transmitted inten- 
sity. The Rayleigh distribution describes large fluctuations in intensity at the output of a 
scattering medium excited by monochromatic radiation. The intensity of a single polariza- 
tion component normalized to its ensemble average has a negative exponential distribution, 
exp(— I / (I)) []J. This distribution obtains under the assumption that the field E = Ae 1 ^ 
can be represented as a superposition of uncorrelated partial waves. However, in multiply- 
scattering media, the coherent nature of wave propagation inevitably leads to both short 
and long-range intensity correlation These give rise to enhanced fluctuations in the 

intensity ||, total transmission ||-10|, and electronic conductance [ TT|— 13 1 , which have been 
studied intensively in the last decade. The degree of nonlocal intensity correlation is a 
measure of the closeness to the localization threshold PHT3|| and determines the statistical 
distributions of key transmission quantities. In this paper we consider fluctuations in pulse 
propagation. This is in contrast to previous studies of the statistics of steady state prop- 



agation in random media and to measurements of the time of flight distribution [14-16 



Fluctuations in the pulse evolution are averaged over an ensemble of samples and one ob- 
tains the arrival time distribution for transmitted photons, which is proportional to the 
path length distribution P(s). This corresponds to the particle transport picture and gives 
a mix of ballistic and diffusive components. Here we consider fluctuations in the dynamics 
of transmission for a given incident and outgoing channel for different realizations of a ran- 
dom medium. Specifically, we define the local or single channel delay time, r ab , as the first 
temporal moment for a transmitted pulse associated with an incident pulse of bandwidth 
Au centered at t = in the time domain and at ujq in the frequency domain, 

J\E ab (t-u ,Au;)\ 2 tdt 

T ab {LV , Au) = . , 1 

J \E ab (t;uJo, ALu)\ 2 dt 

where E ab (t; c^o, Au;) is the transmitted field in channel b arising from an incident wave in 
channel a. This definition has been used in an earlier work in the context of nuclear physics 



17]. This is in contrast to previous discussions of the Wigner time delay and the Wigner- 



Smith time delay matrix, which have been powerful concepts for a statistical description 
of scattering [p!8|-p0|. The diagonal elements Q aa of the lifetime matrix Q = — iS^dS/du 
where S is the 2N x 2N scattering matrix, are interpreted in terms of the time spent in the 
scattering region by a wave packet incident in one channel. As shown by Smith they are 
the sum over all ouput channels (both in reflection and transmission) of the local delay time 
1 2~2j weighted by the probability of emerging from that channel I ab /2N: Q aa = ^ J2 b I a b r ab- 
The sum of the Q aa over all 2N channels is the Wigner time delay t w = J2 a Q aa which is 
the trace of the lifetime matrix and is proportional to the density of states p4|-EH] . Local 
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delay times have been considered for electrons tunneling through barriers and for classical 
evanescent waves j^^^HJ- In these cases, scattering through a fixed structure into a one 
dimensional system is considered. Here we study the propagation of an incident spatial 
mode into a multichannel random medium. The field is detected at a point in the output 
speckle pattern for an ensemble of random configurations. The subscript ab which indicate 
the input and output channels is omitted in the following to simplify notation. 

In a homogeneous medium of thickness L with phase velocity v, the phase accumulated 
as the angular frequency is increased by Auj is A<p = AkL = AujL/u, and is proportional 
to the delay time r = L/v so that the delay time is r = Aip/Auj. In this simple situation, 
the phase derivative 0' = d^/dcu is known as the group delay and is a measure of the 
transit time through a homogeneous medium. In a random medium, however, the wave is 
multiply scattered and the output field is the superposition of partial waves arriving at a 
point. If we ignore fluctuations in the phase velocity for different paths in the medium and 
additional phase shifts associated with focal points and reflections, the phase accumulated 
by the wave in following a path of length s as the frequency is incremented by Auj is 
A<p s = Aus/v, giving a delay time A<j) s /Au> = s/v. The sum of these times weighted by 
P(s), J P(s)s/vds, only gives the average value of the local delay time, (r). In this paper, 
we investigate the fluctuations of r. Measurements of microwave radiations propagating 
through random samples show that the single channel delay time varies with the width of 
the incident pulse. We find that fluctuations in r increase as the pulse bandwidth narrows. 
For narrow bandwidth pulses, large positive and negative values of r are found. These large 
values are related to the nature of the speckle pattern at the output of the medium. They 
occur most commonly when a null in the speckle pattern passes near the detector as the 
frequency is varied. The phase is undetermined at these nulls and jumps by rr radians when a 
null passes through the detector as the frequency is tuned. In contrast to the characteristics 
of static speckle patterns investigated by Freund which are independent of the scattering 
medium | 29fl , here the phase variation reflects the underlying dynamics. When considering 
pulse propagation, it is natural to define the energy transmission coefficient e a b as well as 
the local delay time r a b of the ouput pulse. This gives a set of variables (e, r). We will study 
pulses of specific width, particularly in the limit, Auj — > 0. In this limit corresponding to 
long pulses, the local delay time r approaches the spectral phase derivative 4>' []30| and the 
energy transmission coefficient e approaches the transmission coefficient I. 

The sample studied is composed of randomly positioned |-inch polystyrene spheres at a 
volume filling fraction of 0.52 contained within a one meter long, 7.6 cm diameter copper 
tube. New sample configurations are created by rotating the tube about its axis. Wire 
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antennas are used as the emitter and detector at the input and output surface of the sample. 
A Hewlett Packard 8722C vector network analyzer performs a measurement of the microwave 
field, giving its amplitude and phase. Measurements are made between 7 and 25 GHz, using 
frequency intervals of 625 kHz. This wide frequency range allows us to reconstruct via 
Fourier transformation the time response to an input pulse of any given shape over a wide 
frequency range. The value of the local delay time depends upon the spectrum of the incident 
pulse. Here we consider incident pulses which are gaussian or rectangular in the frequency 
domain. The gaussian pulse has a rapid fall-off in both the time and frequency domains, 
whereas the rectangular pulse allows us to select a precise spectral range but oscillates in 
the time domain. Direct dynamical microwave measurements are possible in principle but 
generating precisely shaped pulses is not always practical. The complex response to an 
incident pulse with carrier frequency u can be written as E(t) = \E(t)\exp(u)ot + <f)(t)). An 
example of the amplitude \E(t)\ and phase <f>(t) of the response to a pulse constructed by 
Fourier transforming the field spectrum in a particular sample configuration is presented in 
Fig. [I] for two pulses centered at 10 GHz with a gaussian envelope l/v / 27rcrexp(— t 2 /2a 2 ) 
where a = 1 ns and 100 ns. 

Various properties have been used to characterize the travel time of a wavepacket ]31 



When the pulse is not appreciably distorted in transmission, the delay time of salient features 
such as the peak can be used. In multiply-scattering media, however, the shape is generally 
unrelated to that of the incident pulse and changes with configuration (see Fig. |I|a). It is 
therefore not possible to associate features of the transmitted pulse with the incident pulse. 
However, the shift of the barycenter of the transmitted pulse intensity at the output surface 
of the sample, as given is Eq. ([]]), is well defined even in a multiply-scattering system. 
We note that the integration over the pulse is reminiscent of the definition of the associated 
quantity, the energy transmission coefficient e(uj , Au) = f \E(t; u , Au)\ 2 dt, where the time 
origin is taken at the center of the incident pulse at the input surface. It is therefore natural 
to use Eq. ([J) to represent the dynamic fluctuations of pulses in mesoscopic systems. We 
find t = 41.8 ns and 59.8 ns for the pulses of Fig. |l|a and Fig. |I]b, respectively. The local 
delay time averaged over 581 sample configurations is (r) = 42.6 ns for a 1 ns incident pulse, 
while (r) = 45 ns for a 100 ns incident pulse. The difference in averaged values is due to 
variations in dynamical properties over the bandwidth. For comparison, the travel time 
through one meter of free space would be 3.3 ns. 

By changing the central frequency of the incident pulse, one can follow the variation with 
frequency of the single channel delay time for a given sample configuration and a given pulse 
bandwidth. The frequency dependence of t(uj, Au) for a square pulse is plotted in Fig. || 
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between 11 and 12 GHz for three different values of the incident pulse bandwidth Aoo. A 
comparison with 

A0. . <f)(uj + Aoo/2) -Mu- Auo/2) 
A^ H = Au; ' 

shows that r and Aifi/ Au do not coincide as expected in such a medium. However, when the 
average of these quantities over sample configurations is taken, (r) is indistinguishable from 
(A<f)/Auj) which is shown as the thin solid line. This figure also shows that fluctuations of 
t and A(p/Auj around their average values are of the same order of magnitude and increase 
with decreasing pulse bandwidth. The local delay time r can be smaller than the travel time 
in free space and even negative. This occurs most frequently for pulses of duration greater 
than (r) corresponding to Au; less than the field correlation frequency, 8u ~ Fig- H 

shows the response to a gaussian pulse with o = 100 ns, centered at 11.1025 GHz, which 
corresponds to the first negative peak in Fig. 0c. 

In order to clarify the character of these fluctuations, we plot in Fig. || the phase derivative 
of the transmitted field, its phase modulus 2tc, and the logarithm of the transmitted intensity 
which represents the complete field between 10 and 10.5 GHz. Large positive and negative 
peaks in the phase derivative are associated with small values of intensity. Indeed a zero in 
the amplitude would correspond to an undefined phase since the real and imaginary part 
of the field are then zero. The equiphase line map around such a phase singularity is a 
"star" and the phase circulation around this singularity is an integer multiple of 2ir 



| 32fl . A phase singularity is by convention positive if the phase circulates counterclockwise. 
Connection between equiphase lines of different phase singularities have been investigated 
by Freund [2Sfl. To explore excursions in the phase as a singularity moves near a point, we 
measured the phase at closely spaced points along a line as the frequency is tuned. After 
the spectrum at a given position is taken, the detector is translated by Ax = 1 mm on a 4 
cm-length line running symmetrically about the center of the output surface. The increment 
in phase along the line from 18 GHz is obtained by unwrapping the phase modulus 2n ||33|| , 
which is shown in Fig. |5|a. As the speckle pattern changes with increasing frequency, phase 
singularities may move across the detection line resulting in a 2n phase difference between 
consecutive detector positions giving a +ir jump on one side of the singularity and — n jump 
on the other side. The phase difference between two consecutive positions of the detector is 
presented in Fig. || A 2tt step is the signature of a phase singularity moving between these 
two positions. As Ax goes to zero, this phase difference plot would become a series of sharp 
steps and flat plateaus. The presence of phase singularities results in large fluctuations of the 
phase derivative 0' as a function of frequency and of detector position as shown in Fig. [|b. 
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Strong fluctuations of Acf)/ Auj will occur when Au; is small since it is then not so different 
from (p' (Fig. 0c). For larger Au; (Fig. |2]a), averaging of random variations in phase reduces 
these fluctuations. 

The magnitude of fluctuations of A<p/ 'Auj and r are comparable as are their correlation 
frequencies. In the limit Au; —>■ 0, these become identical. In Appendix A, we show that the 
delay time can be expressed in terms of the intensity and the phase derivative as, 

where ^^(w) is the spectrum of the incident pulse. For one particular sample configu- 
ration, the rhs of Eq. @ is shown in Fig. [7| for a gaussian incident pulse with bandwidth 
Auj = \j1-no = 79.6 MHz and central frequency Uq between 7 and 25 GHz and is indistin- 
guishable from the delay time r of a gaussian pulse with a = 2 ns which is also plotted in 
Fig. |7|. Eq. (0) is found to be accurate experimentally to within 0.1 % (1%) for a gaussian 
pulse with a = 1 ns (er = 100 ns). Equation @ provides a useful shortcut to the computa- 
tion of the local delay time as compared to computing the temporal integral of Eq. (p]) from 
the Fourier transform of the field. 

For pulse bandwidths much smaller than the correlation frequency, which is essentially 
the inverse of the average delay time, the transmission coefficient I is roughly constant over 
the pulse bandwidth. Therefore, r ~ A(j)/Au for a narrow rectangular pulse in the frequency 
domain (Fig. |2]c). Only at frequencies at which the intensity drops rapidly and cannot be 
taken as constant over the bandwidth are fluctuations in A<J)/Auj larger than fluctuations in 
r. In these phase singularity moves near the detector giving a change in phase of 

the order of 7r over a frequency change Au; that can be arbitrarily small and in particular 
smaller than 8u. The shape of a pulse, which encompasses a bandwidth over which the 
intensity changes appreciably, will generally be different from that of the incident pulse as 
is seen in Fig. [| 

Taking the limit Au; — > in both Eqs. (|A4p and (lABl) of Appendix A, one finds, 



lim t(uj , Auj) = (f)'(uj ) . (3) 

Au- »0 

To illustrate this result, the time delay between 10 and 11 GHz is presented in Fig. ^] for 
various pulse widths and is compared to the phase derivative for the same frequency 
range (Fig. |8|e). We find experimentally that for a pulse with standard deviation a = 400 
ns, which corresponds to a bandwidth of Au; = \j1no ~ 0.4 MHz, the time delay and 
the phase derivative are almost indistinguishable. Finally, considering the total energy of 
the output pulse e(w , Aw) = J \E(t\ uj , Auj)\ 2 dt, we find from Parseval's theorem |J that 
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lim e = I. We are able to define a general statistical set of variables in the time domain 

Alu->0 

(e, r) which capture the statistics of dynamics in mesoscopic systems and which approach 
the variables (J, <$') in the limit of narrow bandwidth pulses. 

In conclusion, we have investigated the dynamics of wave propagation through random 
media by considering the energy transmission coefficient and the local delay time (e, r), 
which, in the limit of long bandwidth limited pulses, approach the transmitted intensity and 
phase derivative, respectively (I,(f>'). We demonstrate that r is the integral over frequency 
of Icj)' weighted by the spectral density of the incident pulse. Fluctuations in r are found 
to be particularly large, for narrow bandwidth pulses, when e is small. This generally 
occurs when a phase singularity in the static speckle pattern passes near the detector and 
indicates that the statistics of r and / are related. Measurements and calculations 
of key dynamical distributions and correlation functions, which will be presented elsewhere, 
confirm the interplay between / and <$' . Here we have dealt with single channel quantities. 
When considering the spatially averaged delay time, however, it is appropriate to weight the 
local delay time by the energy transmission coefficient, which is 1(f)' in the limit of narrow 
bandwidth pulses. The sum of 1(f)' over all incident and outgoing channels is proportional to 
the density of states and is in many respects analogous to the conductance, which is the sum 
of / over all input and output channels. We expect that large fluctuations in the spatially 
averaged delay time will occur as a result of spatial correlation in 1(f)' just as enhanced 
conductance and transmission fluctuations arise as a consequence of spatial correlation of /. 
Thus the observations in this paper should form the basis for treating the dynamical aspects 
of mesoscopic physics. 



ACKNOWLEDGMENTS 



We thank Bart van Tiggelen, Fabrice Mortessagne and Marin Stoytchev for useful dis- 
cussions. We are indebted to Narciso Garcia for valuable suggestions, for support and 
encouragement. This work was supported by the Groupement de Recherches POAN 1180, 
the National Science Foundation under Grant Nos. DMR 9632789 and INT9512975, and a 
PSC-CUNY grant. 



7 



APPENDIX A: 

To prove Eq. (|2|), we first express the transmitted intensity as 

\E(t;u Q , Au)\ 2 = / / dutduiE*^ - u , Au)E{uj 2 -u q , Aw) e i(< ^ 2) *, (Al) 



where E(u — ujq, Auj) = g WO) ^(uj)E(oS) is the Fourier transform of the transmitted field for 
a given pulse shape g^^uiiuj). Taking u)\ = uj and uj 2 = uj + rj, we calculate the numerator 
of r in Eq. (0): 



|-E(t;w , Aa;)| 2 t dt 
= J dt J dr) J dute~ ivt E%u-u ,Auj)E(uj + r}-ujo,Auj) 

= J dt j d V j du[i^(e- ivt )]E*(u-oj ,Auj)E(oj + ri-oj ,Aoj). (A2) 
Assuming that g^ 0} ^{ui) is vanishing at infinity, we find after integrating by parts, 



\E(t;u , Au)\ 2 t dt 

= -% J dr) J dte~ int j duE*{u - u , Au)^-E{uj + rj - u , Au) 

f f 9 ~ 

= —i dr]27!-5(r]) / du E (uj — uj , Auj) — E(u> + 77 — uj , Auj) 

duj E*(u -u ,Alu)— E(uj + ri - u ,Auj) . (A3) 
dui 

Writing E = \E\e %< ^^ we obtain 

j \E(t;uj , Auj)\ 2 tdt = 2it J \E(u - uj , Au)\ 2 (j)'{u) du . (A4) 
The calculation of the denominator of Eq. (|l|) is straightforward: 

j \E(t;uj ,Auj)\ 2 dt = 2iT J \E(u - u , Aco)\ 2 duj , (A5) 

which gives Eq. (0). 



S 



REFERENCES 

[1] J. W. Goodman, Statistical Optics (John Wiley, 1985). 

[2] M. J. Stephen and G. Cwilich, Phys. Rev. Lett. 59, 285 (1987). 

[3] S. Feng, C. Kane, P. A. Lee and A. D. Stone, Phys. Rev. Lett. 61, 834 (1988). 

[4] N. Garcia and A. Z. Genack, Phys. Rev. Lett. 63, 1678 (1989); A. Z. Genack, N. Garcia 
and W. Polkosnik, Phys. Rev. Lett. 65, 2129 (1990). 

[5] A. Z. Genack and N. Garcia, Europhys. Lett. 21, 753 (1993). 

[6] M. P. van Albada, J. F. de Boer and A. Lagendijk, Phys. Rev. Lett. 64, 2787 (1990); 
J. F. de Boer, M. C. W. van Rossum, M. P. van Albada, Th. M. Nieuwenhuizen and A. 
Lagendijk, Phys. Rev. Lett. 73, 2567 (1994). 

[7] E. Kogan, M. Kaveh, R. Baumgartner and R. Berkovits, Phys. Rev. B 48, 9404 (1993); 
E. Kogan and M. Kaveh, Phys. Rev. B 52, R3813 (1995). 

[8] Th. M. Nieuwenhuizen and M. C. W. van Rossum, Phys. Rev. Lett. 74 2674 (1995). 

[9] S. A. van Langen, P. W. Brouwer and C. W. J. Beenakker, Phys. Rev. E 53, 1344 
(1996); P. W. Brouwer, Phys. Rev. B 57, 10526 (1998). 

[10] M. Stoytchev and A. Z. Genack, Phys. Rev. Lett. 79, 309 (1997); M. Stoytchev and A. 
Z. Genack, preprint (Cond.matt. 9805233). 

[11] R. A. Webb, S. Washburn, C. P. Umbach and R. B. Laibowitz, Phys. Rev. Lett. 54, 
2696 (1985). 

[12] P. A. Lee and A. D. Stone, Phys. Rev. Lett. 55, 1622 (1985). 

[13] B. A. Altshuler and D. E. Khmel'mtskii, Pis'ma Zh. Eksp. Teor. Fiz. 42, 291 (1985) 
[JETP Lett. 61, 359 (1985)]; B. L. Altshuler, V. E. Kravtsov and I. V. Lerner, in 
Mesoscopic Phenomena in Solids, edited by B. L. Altshuler, P. A. Lee and R. A. Webb 
(Noth-Holland, Amsterdam, 1991). 

[14] J. H. Watson Jr., P. A. Fleury and F. M. McCall, Phys. Rev. Lett. 58, 945 (1987). 

[15] A. Z. Genack and J. M. Drake, Europhys. Lett. 11, 331 (1990). 

[16] A. Ya. Polischuk, J. Dolne, F. Liu and R. R. Alfano, Opt. Lett. 22, 430 (1997). 

[17] M. Bauer, P. A. Mello and K. W. Mc Voy, Z. Physik A 293, 151 (1979). 

9 



[18] E. Wigner, Phys. Rev. 98, 145 (1955). 

[19] V. A. Gopar, P. A. Mello and M. Biittiker, Phys. Rev. Lett. 77, 3005 (1996). 

[20] P. W. Brouwer, K. M. Frahm and C. W. J. Beenakker, Phys. Rev. Lett. 78, 4737 (1997). 

[21] F. T. Smith, Phys. Rev. 118, 349 (1960). 



[22] Following the thesis of L. Eisenbud, the local delay time is defined in [21] as the ap 



pearance of the peak outgoing signal in channel b, after the injection of a wavepacket 
in channel a. 

[23] E. Doron and U. Smilansky, Phys. Rev. Lett. 68 1255 (1992). 
[24] G. Iannaccone, Phys. Rev B 51, 4727 (1995). 

[25] V. Gasparian, T. Christen and M. Biittiker, Phys. Rev. A 54, 4022 (1996). 

[26] For a complete and insightful review see Y. V. Fyodorv and H.-J. Sommers, J. Math. 
Phys. 38 1918 (1997). 

[27] R.Landauer and Th. Martin, Rev. Mod. Phys. 66, 217 (1994). 

[28] E. H. Hauge and J. A. St0veng, Rev. Mod. Phys., 61, 917 (1989). 

[29] I. Freund, N. Shvartman and V. Freilikher, Opt. Commun. 101, 247 (1993); I. Freund 
and N. Shvartman, Phys. Rev. A 50, 5164 (1994). 

[30] P. Sebbah, O. Legrand and A.Z. Genack, in: OSA TOPS on Advances in Optical Imag- 
ing and Photon Migration ed. R. R. Alfano and J. G. Fujimoto, Vol. 2, p. 386 (1996). 

[31] G. Iannaccone and B. Pelligrini, Phys. Rev. B 49, 16548 (1994-1). 

[32] J. F. Nye and M. V. Berry, Proc. R. Soc. London Ser. A 336, 165 (1974); M. V. Berry, 
J. Phys. A 11, 27 (1978); M. V. Berry, in Physics of Defects, edited By R. Balian, M. 
Kleman and J. P. Poirier (North-Holland, Amsterdam, 1981), p. 141. 

[33] P. Sebbah, O. Legrand, B. A. van Tiggelen and A. Z. Genack, Phys.Rev. E 56, 3619 
(1997). 

[34] A. Z. Genack, P. Sebbah, M. Stoytchev and B. A. van Tiggelen, unpublished. 
[35] B. A. van Tiggelen, P. Sebbah, M. Stoytchev and A. Z. Genack, unpublished. 



10 



FIGURES 



FIG. 1. Amplitude of the time response (solid line) to a gaussian pulse (dashed line) with (a) 
a = 1 ns and (b) a = 100 ns, centered around 10 GHz. Total energy of the input pulse, J \E\ 2 dt 
is normalized to unity. The solid curve is the actual scale, while the dashed curve is rescaled. The 
slowly varying component of the phase <fi(t) is shown as the dotted curve. 

FIG. 2. Local delay time (thick solid line) for square bandwidth pulse with bandwidths of (a) 
Auj = 500MHz, (b) Au = 50 MHz and (c) Au = 5 MHz compared to Acp/Auj (dotted line) versus 
frequency. Averaged local delay time over 581 sample configurations is shown in the thin solid line. 

FIG. 3. Amplitude time response (solid line) to a gaussian pulse (dashed line) with a = 100 ns, 
centered around 11.1025 GHz. Solid curve is the actual scale while dashed the curve is rescaled. 
The slowly varying component of the phase 4>(t) is shown as the dotted curve. 

FIG. 4. (a) Phase derivative, (b) phase modulus 2ir and (c) logarithm of the magnitude between 
10 and 10.5 GHz. 

FIG. 5. Surface plot of (a) cumulative phase and (b) phase derivative with frequency between 
18 and 18.1 GHz across the output of the sample obtained by sweeping the detector over 40 mm, 
each millimeter. 

FIG. 6. Phase difference between 4 consecutive detector positions (x = 19 to 22 mm) at the 
output surface of the sample of Fig. || in the frequency range 18 to 19 GHz. 
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FIG. 7. Right hand side of Eq. (^) between 7 and 25 GHz for one sample configuration when 
considering a gaussian pulse shape of 79.6 MHz bandwidth. This plot is indistinguishable from the 
plot of the delay time versus frequency of a gaussian pulse with a = 2 ns. 

FIG. 8. Time delay between 10 and 11 GHz for an input pulse with (a) a = 10 ns, (b) a = 30 
ns, (c) a = 100 ns, (d) a = 300 ns. The phase derivative for the same frequency range is shown in 
(e). 
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